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Quantum Graphity (QG) is a model of emergent geometry in which space is represented by 
a dynamical graph. The graph evolves under the action of a Hamiltonian from a high-energy 
pre-geometric state to a low-energy state in which geometry emerges as a coarse-grained effective 
property of space. Here we show the results of numerical modelling of the evolution of the QG 
Hamiltonian, a process we term “ripening” by analogy with crystallographic growth. We find that 
the model as originally presented favours a graph composed of small disjoint subgraphs. Such 
a disconnected space is a poor representation of our universe. A new term is introduced to the 
original QG Hamiltonian, which we call the hypervalence term. It is shown that the inclusion of a 
hypervalence term causes a connected lattice-like graph to be favoured over small isolated subgraphs. 

PACS numbers: 04.60.Pp, 04.60.-m 


I. INTRODUCTION 

The idea that space-time geometry may not be a fun¬ 
damental feature of the universe has arisen in several 
speculative models for the microscopic structure of space 
[T], in particular those which aim to formulate a back¬ 
ground independent theory of quantum gravity [2]. This 
approach is taken in research paradigms such as causal 
dynamical triangulation [3], group field theory [1] and 
structurally dynamic cellular networks [S]. It has been 
suggested that the AdS/CFT correspondence in string 
theory also indicates that space-time is not fundamen¬ 
tal [^. To consider geometry to be merely an effective, 
coarse-grained property of space-time may resolve some 
of the conflicts between quantum theory and general rel¬ 
ativity [7]. However, this comes at great technical and 
conceptual cost. A quantum theory with no inherent no¬ 
tion of geometry cannot be formulated with respect to a 
fixed manifold, nor can it make use of the usual global 
symmetries. For this reason, models of emergent geome¬ 
try often use the tools of combinatorics. 

One such model is Quantum Graphity (QG) [5], which 
postulates that space is represented as a dynamical 
graph, where points in space are represented by vertices 
of the graph and adjacencies in space are represented by 
edges in the graph. Since QG uses simple, undirected, 
unembedded graphs, there is no inherent notion of ge¬ 
ometry at the fundamental level in this picture. Rather, 
it is postulated that this model will exhibit two ther¬ 
modynamic phases: a high temperature phase in which 
there is no sensible notion of geometry, and a low tem¬ 
perature phase in which geometry emerges as an effec¬ 
tive phenomenon. The very early universe is believed to 
have been in the high temperature phase, undergoing a 
phase transition at some later time to the low tempera¬ 
ture phase. This process is called geometrogenesis. Since 
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we are interested in examining an emergent property of 
the low-temperature condensed phase, it is natural to 
draw parallels with condensed matter physics, and apply 
intuitions developed there to the problem at hand. In 
this way, QG can be viewed as a member of the newly 
developing family of condensed matter analogue mod¬ 
els of spacetime, which treat space or spacetime itself as 
though it were a many-body system [S]. 

The interpretation of temperature in QG is not clear. 
If geometry is considered to be a collective phenomenon, 
and we are to appeal to statistical physics to explain it, 
then the role of temperature will be significant. How¬ 
ever, when the system under investigation is space-time 
itself, temperature is difficult to define. In this work we 
interpret matter and energy fields living in space-time 
as a heat bath. The graph can lower its energy by cre¬ 
ating matter. This process is not considered in detail 
here, however the string-net condensation of Levin and 
Wen [101 111] has been proposed as a mechanism by which 
matter may be emergent from a QG space-time [Sj. 

The low temperature phase of QG is assumed to be 
a lattice graph, with the geometry of the lattice corre¬ 
sponding to the emergent geometry of space-time [8]. 
In particular, parameters of the model are often cho¬ 
sen so that a regular 2-dimensional honeycomb graph is 
favoured as the ground state. This is clearly not a repre¬ 
sentation of our 3-dimensional spatial geometry, but is a 
useful simplification for studying some general features of 
the model. The high-energy pre-geometric state has usu¬ 
ally been considered to be the complete graph, in which 
every vertex is connected to every other. This graph has 
no geometrical interpretation, making it a candidate for 
pre-geometry. Furthermore, thermodynamic arguments 
support the notion that higher temperature graphs have 
more edges, making the graph with the most edges pos¬ 
sible a natural maximum temperature limit [13] . How¬ 
ever, the empty graph, in which there exist vertices but 
no edges, so that every point in space is completely iso¬ 
lated, is also a candidate high-temperature pre-geometric 
graph [14] . This is perhaps more in line with the concept 



2 


of a universe arising spontaneously from “nothing”, and 
offers a stronger concept of “spacelessness”. 

The geometrogenic phase transition in this model is 
of great significance and has been investigated from a 
number of perspectives. When recast on a line graph, 
the model can be reduced to an analogue of an Ising 
model, making it amenable to a mean-field treatment. 
In this approach, using the average degree of the graph 
as an order parameter, it was found that as the number of 
vertices goes to infinity, the critical temperature of the 
geometrogenic phase transition goes to zero [T3]. An¬ 
other approach has been to estimate the partition func¬ 
tion of the model by considering which types of graphs 
will contribute most m- The two main competing fam¬ 
ilies of graphs are homogeneous lattice-like graphs (due 
to their low energy and consequent high Boltzmann fac¬ 
tor) and random graphs (despite being higher in energy, 
the large number of random graphs means that they 
contribute significantly). Phase transitions between the 
two states (lattice-dominated and random-dominated) 
are found once again to occur at zero temperature in the 
N ^ oo limit. This behaviour is also seen in other related 
graph models which exhibit geometrogenesis min]- A 
zero-temperature phase transition may imply that the 
phase transition does not exist, indicating that the model 
is faulty. Alternatively, it could mean that geometroge¬ 
nesis is a quantum phase transition, rather than a ther¬ 
mal phase transition, and must be modelled accordingly. 
Since its description in QG, the event of geometrogen¬ 
esis has been shown to have a concrete realisation in 
causal dynamical triangulations, and may also be present 
in loop quantum cosmology m- 

In this work we study the growth (ripening) of grains 
from the empty graph under the epitaxial approximation 
(discussed below). We have assumed that the formation 
of spatial grains from a pre-geometric graph follows a 
process analogous to the growth of crystals. Instead of 
a homogeneous transition across the entire graph, we ar¬ 
gue that as the graph lowers its energy it will tend to 
form distinct grains which then knit together to form a 
geometric space. Over time the size of these grains may 
increase until they are sufficient to house the Universe. 
This picture for the evolution of the graph lends itself to 
the formation of domain boundary defects, which have 
been shown to have observable consequences m- 

Numerical simulations of ripening under the epitaxial 
approximation reveal a tendency for the model to form a 
graph consisting of small isolated subgraphs in the low- 
energy regime. Evolution from the empty graph under 
the epitaxial approximation leads to a disjoint graph as 
the lowest energy state of the system, and we show that 
such disjoint graphs are inevitably lower in energy than 
connected graphs, given the Hamiltonian and parame¬ 
ters first proposed for this model [5]. The formation of 
a disjoint low-energy graph was first noted in [14], and 
here we demonstrate that it remains even under more 
accurate calculations. A disjoint ground state is clearly 
undesirable for a model that aims to reproduce an ex¬ 


tended geometric space. Therefore, we propose a modifi¬ 
cation to the model which may lead to the formation of 
a connected space. 

II. MODEL 

We begin by postulating that space may be represented 
by a simple, undirected graph Q, which is defined by two 
sets V and E. V is the set of vertices (or nodes) on the 
graph, and E is the set of edges, which are unordered 
pairs of the elements of V. Two vertices in V, Vi and Vj, 
are said to be connected (or adjacent) if there exists in E 
an edge This graph is dynamic, in the sense that 

the set of edges may change in time, however the set of 
vertices is taken to be fixed. 

The energetics of the graph are determined by the 
Hamiltonian, 

H = Hv + Hl + H^op- ( 1 ) 

The valence term Hy assigns an energy to the graph 
based on the valence, or degree, of each vertex, 

= ( 2 ) 

i 

where Vi is the valence of the vertex, gy is a positive 
coupling constant, p is a dimensionless real number and 
vq is the “ideal” valence of the graph. Hy is minimised 
when every vertex in Q has Vi = vq. 

The loops term Hl reduces the energy of the graph 
when there are more loops, 

■^max 

= ( 3 ) 

L—3 i 

where P(yi,L) is a function that counts the number of 
loops of length L that pass through vertex Vi, r is a di¬ 
mensionless real number and gr is a positive coupling 
constant. The sum over loop lengths L begins at 3, 
because that is the length of the shortest possible non¬ 
retracting loop. Ideally, the sum would extend upwards 
to include loops of infinite length, but to the make the 
model computationally tractable, loop counting is trun¬ 
cated at some maximum length Umax- The weighting 
factor r^/D. is small both when L is small and when L 
is large. Between these points it reaches a peak at some 
value L*, which is determined by r. Thus arbitrarily long 
loops contribute a negligible amount of energy, justifying 
the use of a truncation length Umax, and by varying r we 
can tune the Hamiltonian so that loops of some desired 
length L* contribute most. As r is varied, we find that 
5-loops are dominant over 6-loops for values of r less than 
6, and that 7-loops are dominant over 6-loops for values 
of r greater than 7. Between these points, L* = 6, so 
that 6-loops contribute most significantly to Hl- There¬ 
fore, we have chosen a value for r in the middle of this 
region, r = 6.5. This differs from the value of r = 7.1 
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which was used in previous literature [H]. Repeating our 
calculations with each of these two values of r shows no 
significant difference to our results except where stated. 

The negative sign in means that this term lowers 
the energy of the graph, so that it favours graphs with 
many loops of length L*. The loop counting function 
P(a, L) is determined so that each unique loop is counted 
only once. If such a procedure were not used, then a 
symmetry factor of 1/(2L) would be required to account 
for non-unique loops. 

Rather than count all loops in the graph, we only con¬ 
sider shortest-path (SP) loops, as these are the most rele¬ 
vant for characterizing a lattice-like graph with emergent 
geometry. SP loops are those that contain no “short¬ 
cuts” , so that the distance between two vertices along the 
loop is equal to their distance on the graph. While there 
exist explicit algebraic formulas for counting the number 
of loops on a graph up to loops of length 7, there exist 
no such formulae for SP loops. SP loops are therefore 
calculated algorithmically using the method presented in 

m- 

Finally, i7hop is a kinetic term that allows edges to 
propagate through the graph, changing the configura¬ 
tion of the graph. Here we are not concerned with the 
particular form of this term, but include it so that the 
configuration of the graph may evolve in time. In gen¬ 
eral, adding, deleting or hopping edges does not conserve 
energy. This can be interpreted as the graph being in 
contact with some heat bath, although when our graph 
represents space itself the idea of an external heat bath is 
problematic m- We follow previous work [T^] and take 
the heat bath to represent the creation and annihilation 
of matter on the graph, so that the lowing of the energy 
of the graph can be interpreted as the creation of mat¬ 
ter and the total energy of the graph -|- matter system 
is conserved. Since we neglect matter in this work, the 
evolution of the graph is a non-unitary process. 

As discussed in the Introduction, we take the empty 
graph as the initial state of the system. 


III. RIPENING UNDER THE EPITAXIAL 
APPROXIMATION 

The epitaxial approximation, introduced in m , as¬ 
sumes that evolution of the graph happens at only one 
edge at a time. First, we take the state of the graph to be 
fixed. Then, we add or delete a single edge (equivalently, 
we delete or add a single “edge-hole”). This new edge 
or hole will form in the position that leads to the lowest 
total energy for that graph. Then, the new edge or hole 
is frozen in place and a new edge or hole is added. This 
approximation is useful, as finding the absolute ground 
state is equivalent to calculating the energy for 
different graph configurations, which is not computation¬ 
ally plausible. Making use of graph isomorphism could 
greatly reduce the number of graphs one needs to check, 
but finding and identifying graph isomorphisms is a prob¬ 


lem for which no P algorithm currently exists [22j . 

In |14j , it was shown that when the model evolves epi¬ 
taxially on 24 vertices from an empty initial graph, it 
achieves a local minimum 3-regular state which consists 
of four disjoint subgraphs, and a lower-energy 4-regular 
graph which consists of only 3 disjoint subgraphs. This 
indicates that the graph evolves from the empty start¬ 
ing point by forming small grains and merging them to¬ 
gether, suggesting a process similar to the ripening that 
can be seen in many familiar material systems. 

Consider a material system which is entirely in one 
phase, a, which is suddenly brought into the region in 
its phase diagram of co-existence of two phases, a and j3 
(as may be achieved by a sudden change in temperature 
or pressure). The material will then be in a metastable 
state. Small regions of phase fi will begin to nucleate and 
form particles. If these particles are fairly diffuse and able 
to move free through phase a, we may see the growth of 
the P phase through a process called Ostwald Ripening 
m- Atoms at the surface of a particle are less tightly 
bound than atoms in the bulk, so they are more eas¬ 
ily detached. Once detached, they diffuse through phase 
a until they re-attach to another particle of phase p. 
Since the surface-to-volume ratio is larger for small par¬ 
ticles, these are more likely to lose atoms and less likely 
to gain them, so this leads to larger particles growing 
while smaller particles disappear. In this way, grains of 
phase P form, until equilibrium between phases a and P 
is achieved (or, if the relevant parameters are changed 
further, grains of P will grow until they constitute the 
entire system). 

The process of ripening in QG is quite different from 
that in material processes. The most striking difference 
is that Ostwald ripening occurs due to atoms being less 
tightly bound at surfaces than in the bulk of a material. 
Graphs in QG are non-oriented, so the usual intuition of 
a “surface” does not exist. However, one may define a 
“surface” in QG as a region of vertices which deviate from 
the ideal valence (or, away from the ground state, have a 
larger than average value of |ui — r^ol). If one accepts this 
definition, then the difference in binding energy between 
vertices in the bulk at vertices at the surface is negligible. 
To be clear, neglecting the loops term in the Hamiltonian, 
the energy required to reduce the valence of a vertex in 
the bulk to zero is 

_ ^p(vi-voP 

9v 

where Vi is the valence of the vertices within the bulk. 
The energy required to reduce the valence of a vertex in 
the surface to zero is 

_ ^pvl _ ^p{\vi-vo\ + lp ^ 

9v 

On the other hand, the energy required to promote a 
vertex from the bulk to the surface of a grain is 

_ ^p(|.u._t,Q| + l)2 _ ^p{vi-VoP ^0^ 

9v 
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To make this concrete, we calculate these quantities us¬ 
ing the parameters of Konopka et al [8] and taking 
the valence of the bulk vertices to be 2. This gives 
AEb = 4.9017xl0^ AEs = 4.8899xlO^ and AEb^s = 
118.1903. Here it can be seen explicitly that the energy 
difference between the bulk and the surface is negligible 
compared with the energy required to reduce the valence 
of either of them to zero. 

The model will grow in way that reduces the number of 
vertices which have the maximum value of — ug |, which 
means that when starting from the empty graph the first 
stage of evolution is to distribute disjoint edges through¬ 
out the graph. In the next stage of evolution some 
of these disjoint edges connect to form small “grains”. 
Eventually we have a graph which consists of many small 
disjoint subgraphs. It is from this stage that the process 
we call ripening can begin. 

Consider the case of a fc-regular graph consisting of sev¬ 
eral disjoint subgraphs where k < vq. There are two pos¬ 
sibilities for adding a single additional edge: the new edge 
will either connect two vertices already within the same 
grain, or it can connect two formerly disjoint grains to 
each other. Which one is favoured generally depends on 
the definition of loops we employ in EIl ■ Loop-counting 
based exclusively on shortest-path loops favours an edge 
which connects two formerly disjoint grains, as an inter¬ 
nal edge within a grain will destroy some loops already 
in place. More general loop counting, in which all closed 
paths contribute to Hb (as was used in [Hj), favours in¬ 
ternal edge formation as this creates more loops, whereas 
connecting two separate grains does not (any walk that 
begins in one grain, crosses to the other and then crosses 
back again to the initial grain must traverse the newly 
formed edge twice, so such a walk cannot be a path and 
therefore cannot contribute even in this more general 
loop counting definition). 

The situation from the complete graph is more diffi¬ 
cult to interpret, as it is not easy to see which loops exist 
on the edges when looking only at the holes. In princi¬ 
ple, evolution from the complete case should be able to 
give rise to grains of geometric space, as well as phase 
separation where some regions are in the pre-geometric 
state while others exhibit local emergent geometry. How¬ 
ever, when the graph is already connected the different 
grains will be difficult to identify. Furthermore, algorith¬ 
mic loop counting is computationally prohibitive when 
dealing with highly connected graphs. 

The evolution from the complete graph to a low energy 
ug-regular graph under the epitaxial approximation was 
simulated for N = 24 and N = 36 using parameters 

from Konopka et al [8] {gv = 1, 'Co = 3 and 

DUU 

p = 1.2), the r = 6.5 as discussed above. Important steps 
in the evolution of the N = 36 are presented in Fig[2 

The general features of the evolution presented in Fig[^ 
are mostly the same as those that were seen in the N = 
24 case. The only significant difference is the fact that 
the ug-regular state on TV = 24 consists entirely of the 
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FIG. 1. a) The epitaxial growth of spatial domains start¬ 
ing with the empty graph as the initial graph, b) The graph 
always evolves in a way that reduces the nnmber of vertices 
with the maximum value of |ui — uo|, so the first step of evo¬ 
lution is the formation of isolated edges, c) Isolated edges 
then form 4-loops, d) which knit together in a process we call 
“ripening” to form cubic structures. The graph in d) is a local 
minimum, and is lower in energy than a connected, isotropic 
honeycomb lattice for the values of the parameters we have 
used. 


symmetrical cubic structures seen in Figj^d), whereas it 
is not possible to populate a graph on 36 vertices entirely 
with cubes as 36 is not divisible by 8, so a small K 4 forms 
from the “leftover” vertices. Thus we have a frustrated 
ground state. Frustration is expected to be negligible in 
the N ^ 00 limit. 

Somewhat surprisingly, the formation of isolated 4- 
loops in step c) of Fig is more favourable than the 
formation of isolated 6-loops. This is because, although 
6-loops reduce the energy more, the graph is able to pro¬ 
duce more 4-loops, so the overall effect is a lower energy. 
Furthermore, epitaxial growth will favour the formation 
of 4-loops as these minimize the instantaneous energy 
at each step in evolution. Even more surprisingly, the 
ground state configuration in Figj^d) consists of isolated 
cubes, and not a connected honeycomb lattice. This cu¬ 
bic state is lower in energy than a connected, isotropic 
honeycomb lattice for the same values of the model pa¬ 
rameters, the honeycomb having an energy of 29.4048 
while the isolated cubes have an energy of 29.0779. This 
difference is small, but it is enough to demonstrate that 
the honeycomb graph is not the ground state, and the 
energy gap is expected to grow with N (see Fig. and 
the surrounding discussion). The honeycomb graph does 
contain more 6-loops than the isolated cubes (the honey¬ 
comb has 2/3 loops per vertex, while the isolated cubes 
have 1 /2 in the ideal case and 2/9 with the K 4 defect seen 
in Fig. [^d), however the cubes also have many 4-loops. 
It is possible to increase the value of the Hamiltonian pa¬ 
rameter r to a value such that the 4-loops become neg¬ 
ligible. For example, the honeycomb is lower in energy 
than the isolated cubes for r = 7.1 (the value used in 
previous work on QG 0), however at this value 7-loops 
dominate over 6-loops. If a connected regular graph were 
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FIG. 2. The evolution of the graph from one consisting of 4- 
cycles seen in Fig[^c) to one constituted by cubic structures 
seen in Fig d) is presented step-by-step. Each subfigure 
in sequence represents the formation of one additional edge. 
Interestingly, the diameter of each grain increases from 2 in 
a) to 4 in b), before decreasing back to 3 in c). This is in 
contrast to the material version of Ostwald ripening, where 
the average diameter of large particles grows monotonically. 


to be favoured at all with this value, it would be a graph 
of heptagons, which would constitute a hyperbolic rather 
than a flat geometry. The crucial consequence of this dis¬ 
cussion is that the connected honeycomb graph cannot 
be the ground state of the model, and this generic effect 
of favouring disconnected over connected graphs is likely 
to hold for other parameter choices. 

The process of ripening happens primarily between 
steps c) and d) in Fig Two “squares” in c) knit 
together to form “cubes”. The cubes form one-by-one 
under the epitaxial approximation. The formation of a 
single cube from two squares is shown in Fig[^ 

In more traditional descriptions of ripening in mate¬ 
rial systems, one of the key quantities of interest is the 
average radius of particles as a function of time (as well 
as other factors such as the coefficient of diffusion). For 
abstract graphs with no inherent notion of geometry, the 
radius of the graph is defined as the smallest eccentricity 
of the graph, the diameter the largest eccentricity, where 
the eccentricity of a vertex is the largest possible geodesic 
distance to any other vertex on the graph. The case of 
ripening in QG as shown in Fig is interesting because 
the diameter of the grains does not increase monotoni¬ 
cally. Rather, we see the diameter of the grains undergo 
a sharp increase from 2 to 5 between steps a) and b), 
followed by a decrease to 3 in step c) (where it remains). 
This “two steps forward one step back” type of growth 


is also seen when the graph transitions from one made 
up of isolated edges to one consisting of isolated 4-loops. 
The radius of the grains, however, does increase mono¬ 
tonically from 2 to 3. If the graph theoretic concepts 
of radius and diameter apply to the emergent geometry, 
this implies a picture where space forms by expansion, 
followed by slight contraction. However, the extent to 
which the graph theoretic notion of diameter of the grains 
corresponds to the geometric concept of the diameter of 
a space is unclear. 

Another way in which the evolution observed here dif¬ 
fers from ripening in material systems is that the process 
occurs mostly homogeneously. The 4-cycles in Fig die) 
form one at a time, but cubic structures do not begin 
to form until there are no remaining vertices of degree 
1. Likewise, the 4-cycles themselves do not begin to 
form until there are no remaining vertices of degree 0. 
The graph goes through distinct stages of evolution, but 
goes through them uniformly. This is in contrast to the 
nucleation and phase separation that is observed when 
material systems undergo ripening. This feature is most 
likely a consequence of both the epitaxial approximation, 
and the form of iLy, which causes evolution to focus on 
vertices of the maximum value of \vi — z)o| above all else. 

Our simulation shows that the formation of two sepa¬ 
rate grains in this model is preferable to a regular lattice 
spreading out from the boundary of a honeycomb grain. 
We may then expect that a QG model with the empty 
graph as a high-energy starting point tends to form dis¬ 
joint subgraphs. Such disconnected states can be lower 
in energy than a connected, extended lattice. Indeed, 
the cubic graph seen in Fig[^d) is lower in energy than 
an isotropic connected honeycomb lattice on the same 
number of vertices. This undermines an assumption of 
the model that the Hamiltonian in equation [^favours an 
extended honeycomb lattice. 

Honeycomb lattice grains have the topology of a flat 
torus, so that there are no “loose ends” (no vertices with 
Vi ^ vq). This can be thought of as a rectangle with 
opposite edges identified, tiled with a honeycomb lattice. 
The “lengths” of the two sides of this rectangle (alterna¬ 
tively, the minor and major circumference of the torus) 
are important, as when the sides are of length 6 then 6- 
loops can be formed by circumnavigating the space, not 
just by the plaquettes in the lattice. This increases the 
number of 6-loops and thereby decreases the energy. So, 
for a grain of a given N which we assume to be con¬ 
nected, the lowest energy honeycomb conhguration will 
be a torus which is long and thin, so that 6-loops can 
wrap around the width of the torus. 

In hgure we show the energy of various graphs as a 
function of N, as calculated with algorithmic loop count¬ 
ing. The blue curve shows graphs which are flat tori with 
a minor circumference of 6, so that 6-loops may form by 
winding around the width of the torus. Therefore N is 
changed by changing the major circumference, i.e. mak¬ 
ing the torus “longer”. The red curve shows the energy 
per vertex of isotropic toroidal graphs, where the ma- 
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FIG. 3. Plot of the energy per vertex of various graphs 
as a function of the number of vertices N, using parameters 
based on Konopka et al [8]. The red line shows the energy 
per vertex of isomorphic toroidal honeycomb graphs, where 
the major and minor circumferences are equal. The blue line 
shows honeycomb graphs where the minor circumference is 
fixed at 6, and only the major circumference is varied. For 
comparison, with our parameters a graph made up of isolated 
cubes always have an energy per vertex of 0.6385, which is 
notably lower than any of these honeycomb graphs. From 
this figure, the N ^ oo limit is not certain, but there does 
not seem to be any indication that larger subgraphs will be 
more favourable at large N. 


jor and minor circumferences are equal. It can be see 
that the 6-by-6 torus does indeed have the lowest energy 
per vertex of the honeycomb graphs shown, with the en¬ 
ergy per vertex of the long thin tori becoming constant 
for large N. A graph made up of isolated cubes has an 
energy per vertex of 0.6385 with our parameters, which 
is significantly lower than any of the honeycomb graphs 
presented in Fig[^ This shows a clear tendency for the 
formation of a disconnected space, which would make 
the model a poor description of our Universe. This is 
discussed further in the Conclusion. 

It is worth noting that, although it is not always pos¬ 
sible to tile an arbitrary space with a honeycomb lattice 
(see, for example, a buckyball, which must include some 
5-loops in order to adopt spherical topology), this is a 
geometric constraint and has no bearing on the topology 
that a graph may adopt. 


IV. HYPERVALENCE 

The existence of a connected lattice-like ground state 
is an important requirement of QG. All of the terms in 
the Hamiltonian and all of the values used in our cal¬ 
culations were originally tailored to favour a connected 
honeycomb graph as the low energy state, however as 


we have seen there exist disjoint states of lower energy. 
The simplest solution to this problem is to modify the 
QG Hamiltonian in Eq by introducing a term which 
will drive connectivity. In particular, here we introduce 
a generalization of the valence term Hy (Eq. , which 
we call the “hypervalence”. 

While the valence depends on the degree of each vertex, 
the hypervalence also depends on the 2-degree and 3- 
degree, i.e. the number of vertices at distance 2 and 
3 from the initial vertex. The hypervalence energy is 
defined 


H 


HV 


= E 


9vi^' 






(7) 

the parameters here are all analogous to those in the 
valence term. is the number of vertices at distance d 
from vertex i, g^d is a positive real number determining 
the strength of the d**' term and pd is a dimensionless 
real number determining the penalty for deviation from 
the ideal d-degree, Vd- 

In principle, the number of terms within the hyperva¬ 
lence term could extend to N. Eq. [^has been restricted 
to three terms for computational tractability, and there¬ 
fore the weighting parameters g^d and pd must be selected 
so as to be negligible for d > 3. Rapidly decaying cou¬ 
plings also implies a kind of pseudo-locality in the model 
- vertices an arbitrary distance away do affect the energy 
contribution from a given vertex, but this contribution 
quickly becomes negligible as the distance increases. 

For the modified Hamiltonian to give rise to a regu¬ 
lar ground state, we require Vd = dvi. In keeping with 
the original model, we set vi = vq = 3, gvi = ffy = 1 
and Pi = p = 1.2. For higher-degree terms to con¬ 
tribute negligibly, coupling constants were chose to scale 
as gvd = and the pd were fixed so that the ar¬ 

gument of the exponential of each term in the hyper¬ 
valence is equal for the empty graph. This leads to 
Ps = ^ Other values are in principle possible, so 

long as higher-degree terms contribute significantly less 
than lower-degree terms. 

Numerical simulation of the ripening of this new model 
was performed starting from an empty graph on Y = 36 
under the epitaxial approximation, with various stages of 
the ripening process along with the final low-energy sate 
shown in Figure]^ This low-energy graph is connected, 
making it a better representation of a continuous space 
than the low-energy graph in figure However it is still 
not the kind of lattice-like graph hoped for by the QG 
program, nor does it have any obvious interpretation in 
terms of an emergent geometry. 

Explicit calculations show, however, that with the 
hypervalence term, an isotropic connected honeycomb 
graph is lower in energy than both the isolated cubes 
and the low-energy state found by the epitaxial approx¬ 
imation shown in Fig. The honeycomb graph has an 
energy of 38.5950, while the graph in Fig.|^f) has an en¬ 
ergy of 46.2645 and the isolated cubes shown in Fig. did) 
have an energy of 5.0830 x 10"^. There are two important 
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FIG. 4. Stages of the ripening process under the hyperva¬ 
lence term. It can be seen that this does favour a connected 
graph, but the low energy state in f) is not lattice-like and 
there is no immediately evident emergent geometry. The low¬ 
est energy state here is still higher in energy than the regular 
2-dimensional honeycomb lattice, indicating a failure of the 
epitaxial approximation to arrive at the ground state of the 
model. 


implications of this: 1) the addition of the hypervalence 
term is sufficient to favour an extended, flat, lattice-like 
ground state, and 2) the epitaxial approximation is not 
sufficient to find the ground state. The insufficiency of 
the simple epitaxial approximation employed here implies 
the importance of higher-order process, where multiple 
edges may be added or deleted simultaneously. However, 
simulating such processes is computationally difficult and 
beyond the scope of this work. 


V. CONCLUSION 

Ripening under the epitaxial approximation leads to 
the formation of a disconnected space. Clearly, a dis¬ 
connected space is not a good representation of the Uni¬ 
verse in which we live. The concept of disconnected space 
may hint at some interpretation in terms of a multiverse 
scheme, but only if any of the individual grains in the 
graph are large enough to support the large extended 
Universe we find ourselves in. This tendency to form dis¬ 
connected spaces was noticed previously m, and was a 
key motivation for this work. 

The epitaxial approximation considered in this work 
only allowed for one-edge processes, where only one edge 


may be created or deleted at any time. Extending this 
approach to include the possibility of two-edge processes, 
where either two edges are created or deleted simultane¬ 
ous, or where two edges hop to different positions simul¬ 
taneously, may radically alter the graph dynamics pre¬ 
sented here. It must be assumed that such processes 
are lower-probability events, but if they are sufficiently 
energetically favourable when compared with single-edge 
processes they may still dominate, or at least influence, 
graph evolution. 

We have demonstrated that the honeycomb graph is 
not the ground state of the original QG model. Since the 
model was designed to give the honeycomb graph as a 
ground state, this may be considered a partial failure of 
the model. However, we have demonstrated that a simple 
extension of the model to include a hypervalence term 
restores the possibility that the connected honeycomb 
graph may be the ground state. Furthermore, we have 
shown that the epitaxial approximation employed in |14j 
is not sufficient to return the true ground state of the 
model, at least when the model is extended to include a 
hypervalence term. 

We have proposed a solution to this connectivity prob¬ 
lem by including an additional term in the Hamiltonian, 
which we have called “hypervalence”. This term depends 
not only on the degree of each vertex but also on the 2- 
and 3-degree (and, ideally, all higher degrees in a rapidly 
converging manner). Hypervalence is a natural extension 
of the valence term in the original model, and its inclu¬ 
sion leads the model to favour a connected graph. While 
isolated cubes were found to have lower energy than a 
connected honeycomb lattice in the original model, the 
honeycomb lattice does become a lower energy state and 
a potential ground state when the hypervalence term is 
considered. However, the lowest energy state found un¬ 
der an epitaxial approach was not the honeycomb graph, 
nor any other lattice-like graph with a clear interpre¬ 
tation as a geometrical space. The lowest-energy state 
attained by this approach was also higher in energy than 
the honeycomb graph. 
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